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Kazufumi Ito, Jeffrey S. Scroggs, and Hien T. Tran 
Center for Research in Scientific Computation f 

North Carolina State University l 

Raleigh NC 27695-8205 , USA ■ . 

The optimal boundary temperature control of the stationary tlermally coupled incompressible Navier- 
Stokes equation is considered. Well-posedness and existence of the optimal control and a necessary op- 
timality condition are obtained. Optimization algorithms based on the augmented Lagrangian method 
with second order update are discussed. A test example motivated by control of transport process in 
the high pressure vapor transport (HPVT) reactor is presented to demonstrate the applicability of our 
theoretical results and proposed algorithm. 

1. Introduction 


In this paper we discuss the optimal control problem of the stationary thermally 
coupled incompressible Navier-Stokes equations. Consider the following optimal 
control problem 

minimize J(g) = <p(u, T - T 0 ) + - \g - T 0 |^ (r|) , g€C (1) 
subject to 

u ■ + V/7 = v An + y (7 - 7b) e d + / , 

V • « = 0, u| r = 0 , 

u ■ VT = V ■ {k V7) , ( } 

T = T 0 on T 0 and n -VT = H (g - T ) on T , , 

where / e L 2 (Q) U is a source field, u, p, T stand for the nondimensionalized 
velocity vector in R d with d = 2, 3, pressure, and temperature, respectively and 
C is the closed convex set in L 2 (ri) such that 

T\ <g <T 2 on r i (3) 

Here, f\ < To < tz and T,, i = 0, 1 are disjoint open sets in T such that T — 

To U Tj . The constant y is given by y = — where To > 0 is a constant reference 

To 

temperature and g is the gravitational constant, and ej denotes the d — th unit vector 
of R d . Throughout this paper we assume that £2 is sufficiently smooth and v, k 
and H are positive constants. 

This control problem is motivated by control of transport and growth processes 
in the high pressure vapor transport (HPVT) reactor [13], For example, we may 
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consider the Scholz geometry (Figure 1 in [13]). The source material and the grow- 
ing crystal are sealed in a fused silica ampoule that is heated by a furnace liner at its 
outer cylindrical surface. The substrate T o (the single crystal) is located on a fused 
silica window (the bottom of the ampoule) which is cooled by a jet of helium gas 
from the outer surface. HPVT processes are based on physical vapor transport and 
can be described very roughly as proceeding via evaporation at the polycrystalline 
source and condensation at the surface of the cooler substrate. The system (l)-(2) 
is called the Boussinesq equations, where we assume tbat the flow is incompress- 
ible and the transport phenomena of a single (carrier) ^s is modeled. At the wall 
we assume Newton’s law of cooling holds. 

The cost-functional can be of tracking type 


i, T - T 0 ) = [ | u- u d \ 2 + \T - T d \ 2 dx 
Jn 


where ( u d , T d ) is the desired state, or minimization of friction force of flow in a 
subregion £2] of £2 

<P(u,T-Tq)= f |V x u\ 2 dx . 

JCi x 

2. Well-posedness 

In this section, we discuss existence and uniqueness of solutions to (2). Let U = 
L 2 (V i) be the control space, V = V 0 x V\ where V 0 is the divergence-free subspace 
of (W 0 l (fi)) a [7] and 

V. = /4 0 = H' (S2): *|r 0 =0} 

and set H = H 0 x L 2 (£2). Ho is the closure of Vo with respect to the L 2 (fi) </ -norm, 
and is defined by 

Ho = {<{> € L 2 {i 2) d : V ■ <p = 0 and n ■ <p = 0 on T). 

H is equipped with the natural L 2 -norm and V is equipped with the norm 


l(0> X)iv = I0lv o + Ixlv, 


where 


Mvb = Wliito, and Ixlv, =|Vxli HQ) + //|xli 2 (ri) 
for ir = ( <p , x) € V. Define the trilinear form b on H x (Sl) d by 

b(4> i. <Pi, <Pj) = (<Pi ■ V (p2 , cfo) 
for <j>i e H l (Q) d , i = 1, 2, 3. Then we have 
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Lemma 1 The trilinear form h satisfies 

(а) |6(0] , ,01, 0 3 )| < |0| ii/* 1^02 It. 2 1 03 It. 4 < M t |0] |//i |02!W |0l|/t < 

(б) bub ] . fiz. 0?) -i- b{<b \ , 03, (pz) = 0 

provided that V ■ 0 t = !) and (n ■ 0 | )(<pz ■ 0 i ) = 0 on V 

(c) 16(01,01,03)1 < /Wli0l ||.:'!0I lJ/r103li^I03li/rl0=l«* for d - 1 

(d) |6(0i, 02 , 0 3 )| < Mn |0i i ^ j0i i ^i^I03 i ^ 4 |03 I0ai w * for d = 3 

/or 0 , € W , (S 2 )‘'. t = 1 . 2,3 
Proof: By Green’s formula 


6(0], 02, 03 } + 6(01, 0 3 , 02) = (0|, V (02 • 0}J) = (« - 01 , 01 • 0 3 )r 

for 02 e C‘(Q) and V ■ 0 ] = 0 . Hence (b) follows from the continuity of 6. The 
last two assertions follow from the fact that \ 5 c |0!/ ; '|0| w r tor d = 2 and 

|0 12,- 1 < c |0|^( 4 |0i^f for 0 6 Z/ 1 (Q) and some constant c. □ 

In particular. Lemma 1 implies that 

6 (u. 0 . 0 ) = O for u . 0 c Vo. ( 5 ) 

The weak or variational form of ( 2 ) is given by 

v (Vu, V0) +b(u. u. 0 ) = (/ -I- y {T - 7 0 )<; u -. 0 ) (6) 

for all 0 6 V 0 and 

K (V(T - 7 0 ). V X ) + (« • V(T - To). Xk-xv, 

+ k (H (T - T 0 ), x)r, =x (H (g - 7q), x)r, • (") 

for all x e V\. The pair (u. T — T 0 ) e V is said to be a weak solution of ( 2 ) if (6), 
( 7 ) holds for all 0 = ( 0 , x) e V. Then we have the following theorem. 

Theorem 2 Given g € L 2 (T 1) there exists a weak solution (u. T — To) e V to 
( 2 ) and 

l(«- T - To)| v < const (|/| L : + Igit^r,)) • 

Moreover, i iff\ < g(x) < fz a.e. in x e T 1 then f\ < T(x) < 7 ? a.e in Q for 
every solution (u.T — To) e V. 

Step I: (Existence) We show that (6), ( 7 ) has a solution ; = (a. T — To) € V. 
Given u e Vo, we consider the linear equation 

v {Vu, V0) -i- 6(«. u, 0 ) — (y (T - To)ej + /, 0 ) for 0 6 V 0 (8) 
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^(V(T-To), V x ) + <5- V(T-To), X> + /c(tf(T-T 0 ),x ) 

= K (ff (* - To), x)r, , (9) 


which implies 


for x 6 Vi . First, we show that (8), (9) has the unique solution («, T — To) e 
V. Then, we show that the solution map S on Vo defined by S(ii) = u where 
(u, T — Tq) e V is the unique solution to (8), (9) has a fixed point by Schauder fixed 
point theorem. The fixed point u € Vo and the corresponding solution T — To € V\ 

define a solution to (6), (7). By Green’s formula ^ 

' « 

(u ■ V$, x) + (u • Vx, 9) = 0 and in particular (u-X>X) = 0 (10) 

for € Vq, and 0, x € //'(ft). Hence, from Lemma 1 the sesquilinear form 

«(Vx,, Vxi) + (u- Vxi, X2 )+(kH xi, X2)r, =0 

on V\ x V] is bounded and V\— coercive. It thus follows from the Lax-Milgram 
theorem that equation (9) has a unique solution T — To € V\. Choosing x = T -To 
in (9), we have (independent of u € Vo) 

|T-T 0 | 2 Vi < H |g - T 0 || 2(ri) . (11) 

Next, the sesquilinear form on V 0 x Vo defined by 

v (V0i , V0 2 ) + b(u, <t >\ , 0 2 ) 

is bounded and Vo— coercive from Lemma 1 and (5). Thus, by the Lax-Milgram 
theorem, equation (8) has a unique solution u e Vo, and we have 

\u\v 0 < — (I/It.* + Y \T - TokO (12) 

v 

where |d>|// 0 < 101 v 0 r <t> € Vo. Let C be a closed convex subspace of V 0 , 

defined by 


and thus \u k - 
there exists at 

Step 2: (L°° e 


for all solutioi 
X € V, [19] a: 

K 

Since from (1 


where 

Thus, we obta 
can prove thai 

We have r 
/ and g - T 0 . 

Theorem 3 1 

unique solutii 

Proof: Suppo 
have 


C = [<t> € V 0 : |0| v 0 < -f (l/ki +yM 4 H |g - Tol^r,)} , 

where |x Ir. 1 £ Ixlv, for x € V\. Then it follows from (11)-(12) that S maps 

from C into C. Moreover, the solution map S is compact. In fact, if u k converges 
weakly to « in V 0 then lu* — u\ L * — ► 0 since H x (£2) is compactly embedded into 
L 4 (£2). Let (ut, T k — To) € V and («, T — 7o) e V be the corresponding solution 
of (8),(9), respectively to u k 6 Vo and a € Vo. Then we have 

* (Tk - T, x)v, + (( u k - 5) • V(T - T 0 ) + u k ■ V(T* - T), x) = 0 


v(V5, 

k(VT 

for 0 e V 0 an 
4> — u and x 


for x € V] from Lemma 1 and (10) 

k | T k — T|v, < M\ \u k — «|i*|T — TqIv, 
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^ax-Milgram 


(12) 


>space of V 0 , 


i that S maps 
Hk converges 
n bedded into 
ding solution 


() =0 


which implies ] 7* — T\v, — > 0. Similarly, we have 

v I u k -u\ Yo < M\ | u k - u\ L *\u\v 0 + yM-i | T k - T\ c . 

and thus|u* — «| Vo — > 0. Now, by the Schauder fixed point theorem (e.g, see [19]) 
there exists at least one solution to (6), (7). 

Step 2: ( L°° estimate) We show that if 7] < g <^F 2 then 

<■ 

T\ < T < 7*2 a.e. x e Q . 

for all solutions (a, T — To) e V to (6), (7). In fact, let x = inf(7\ f\). Then 
X e V\ [19] and we have from (6), (7) 

K (V7\ Vx) + (« ■ vr, x) + (kH (T - g), x)r, = 0 . 

Since from (10) (« • V7\ x) = 0 we have 

ic (Vx, Vx) + (ac// ( T - g), x)r, = 0 

where 

(T ~ g)X = (T — f\ -(g-f,))x >lxl 2 on r , . 

Thus, we obtain |xl \ = 0 which implies x = 0 and hence T >T\. Similarly, one 
can prove that T < to, choosing the test function x = sup(T, to). 

We have also the uniqueness of solutions under the smallness assumption on 
/ and g - T 0 . 

Theorem 3 //|/| t 2 and ~|g — Tol^r,) are sufficiently small then (6), (7) has a 
unique solution in V. 

Proof: Suppose («,-, 7} — To) € V, i = 1, 2 are two solutions to (6), (7). Then we 
have 

v (Vu, Vcp) + («i ■ Vu + ii ■ V«?, <p) = y if ej, i p) 

* (Vf, V X ) + («, • Vf + 5 • V(T 2 - To), x) + Oetf T, x)r, = 0 

for (p e V o and x 6 Vi , where u = «] — h 2 6 V 0 and T = 7) — 7 2 e V) . Setting 
ip = u and x = 7 we obtain from (5) and (10) 

v l«lv 0 < Mi |« 2 l Vbl«l v 0 + KM3M4 |T| V| l«l Vb 

tc (|V7| 2 + // |7| 2 ( ) < M,|T 2 - Tolv, |«k|f|v, • 

If we set X — |m|v„ and Y = \ f\v [ then this implies 

(v - M, \u 2 \ v „)X 2 < yM 3 M 4 XY and k Y 2 < Mi| 7 2 - 7o|y,X Y . 
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Hence if k(v — M\ |u 2 lv 0 ) — yM\ A/3A/4 \T 2 — 7b|v, > 0 then X — Y = 0 and thus 
(«i, T\) = (u 2 , T 2 ). From (11)-(12) we have 

\Ti — 7b|v, < H \ g — 7o| L i (ri ) 

M 3 (13) 

l« 2 lv„ < (1/It, + YM 4 H \g - Toltur,,). 

Thus, if |/| Li and ^|^ — 7bli2(r,) are sufficiently small then (6), (7) has a unique 
solution. O 

£ 

Moreover, we can make the following demonstration of the regularity of the 
solution («, T — To) £ V. Define the Stokes operator A on Hq by 


with domain 


( Au , <t>) Ho = (V«, V0) for 0 e V 0 


(14) 


Then Problem 

Proof: Let ( 

0 > °, I gk 
I («*. T k - Ty 
by the same i 
V x C since 
tri linear forn 


for Zi = (u, 
bedded into 
0 6 V, whe: 
the limit (z, 
lower semici 


dom(A) = {« € Vq : | (V«, V0)| < c|0|// o for all 0 e Vb) . 


(15) 


Then it is known [18] that A is a positive self-adjoint operator on Hq, dom(A) c 
H 2 (£2) 3 and Vq =dom(A ,/2 ) = [Ho,dom(A)]i/ 2 - Let d = 3. Since 

|(m • Vu, 0)1 < |M|t*l VMji.2101/,3 

for u, cp e Vq and 

H' (£2) C L 6 (£2) and V 1/2 C L 3 (Q) , 
where V )/2 = [Vb, Hq] i/2 , uVu € VL 1/2 =dom(A~ 1/4 ) for u e V 0 . Thus, 
u = A _l (p(r-ro)e 3 +/-M-V«) 6 dom(A 3/4 ) = [0, dom(A)] 1/2 C H 3/2 (i 2) 3 . 
Hence, u • Vu e L 2 (£2) d and a edom(A) C // 2 (£2) 3 . 


si 

I 


1 


3. Necessary optimality condition 

We now show the existence of solutions to the optimal control problem (l)-(2). Let 
us denote by S(g), the solution set of (6), (7) for g e L 2 (T 1 ). 

Theorem 4 Consider the minimization problem ( 1 )-(2): 


minimize J(g ) = <p(u, T — 7b) + § \g - T 0 | 2 J(ro 

over ( u , T — To) € S(g) and g e C. 


(16) 


Assume that <p satisfies 

( p{z) : z = (u,T — 7o) € V — »■ (J convex and lower semicontinous 

and <p{z) <b\ |z | \ + b 2 for b\, b 2 e R + . 
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(14) 
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. dom(A) c 

tus, 

H 3 / 2 (Q ) 3 . 

(D-(2). Let 
(16) 

itinous 


Then Problem ( l)-(2) has a solution. 

Proof: Let (( u k , T k — T 0 ),gk ) € 5(g* ) x C be a minimizing sequence. Since 
P > 0. |g;. — Tol^r,) is uniformly bounded in lc and thus from ( 1 1 )-( 1 2) so is 
|(t/*, 7* — T q )\ v . Hence there exists a subsequence of (£}, which will be denoted 
by the same index, such that (a*, 7* —To, gf) converges weakly to (u. T — 7 0 ,g) e 
V x C since V x L 2 (Ti) is a Hilbert space and C is closed and convex. Define the 
trilinear form b on V 3 by 

..«r 

0(Z|, Z 2 , f) = b{U\,U 2 , <j>) + • V02, x) 

for Zi — ( Uj,9i ), i = 1,2, i ]/ = (<p,x) e If. Since //'(£2) is compactly em- 
bedded into L 4 (Q) it follows from Lemma 1 that b(Zk, Zk, ft) -*• b(z, z, i p) for 
f/ e V, where Zk = (w*, Tk — 7o) and z = (u. T — To). Hence, for ft = (0, x) 6 V 
the limit (z, g) satisfies (6), (7) and thus z € 5(g). Now, since is convex and 
lower semicontinuous, it follows from [5] that (z, g) minimizes (16). □ 

Problem (16) is equivalently written as a constrained minimization on x = 
((«, T - T 0 ), g) e X = V x L 2 (T,) with 

minimize J(x) = <p(«, T — 7o) + f |g - 7 0 | 2 J(ri) over.t e X 
subject to c(x) = 0 and g e C 

where the equality constraint e : X — > Y = V* is defined by (6), (7); i.e„ 

(e(-r), y>) = a(z, f) + b(z. z, ft) - k (H(g - 7b), x)r, -(/, 0) (18) 

for i/t = (0, x) € V, where 

a(z. ft) = v (V«, V0)- y ((7— 7 0 ) e d , 0)+* (V(7-7 0 ), V x )+k// (7-7 0 , x)r, 

forz = («, 7 -7o), f/ = (0, x) € V. Assume that x* = (z* = («*, 7*-7 0 ),g*) 
denotes the optimal pair of (16). Then we have 

Theorem 5 Assume that x* is a regular point in the sense [14] that 

0 € int [e'(x*)(u. h — g*) : v 6 V and h e C ] . (19) 

Then there exists a Lagrange multiplier k* e V such that 

a(r , 0) + b(ft, z\ A.*) + b{z‘, ft. A.*) + <*'(z*), 0) =0 (20) 

/or ft e V and 

{Pg m - kHV 2 , h -g*)r, >0 for all h e C. □ (21) 

Proof: It follows from [14] that if (19) is satisfied, then there exists a Lagrange 
multiplier A.* = (A.*, A.£) € V such that 

(<p'(z*), ft) + P (g*, h - g*) r , + e'(x m )(f/, h - g‘) > 0 
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for all it eV and h e C, that is 

(*V). f) + P ig*. h - 8') r, + a(X\ xjt) 

+&W, z\ A.*) + h(z*. xjt, A.*) — tcH (h — g\ X*) r , >0 ' 

for all i/r e V and h e K. Setting xjt = 0, we obtain (21). Next, setting h — g m in 
(22), we obtain (20). O 

Concerning the regular point condition (19), ^e have 

Lemma 6 If g* e int ( C ) then the regular point condition (19) is equivalent to the 
condition that for v = (v \ , 1 ) 2 ) 6 V 

a(v, xjt) + z*. u) + 6(z*, tM) = 0 

for all ft e V and V 2 = 0 on Ti 


implies v = 0. 

/Voo/- If g* e int (C) then (19) is equivalent to the condition that G = e'(x *) is 
surjective. Define the linear map C e £(X, V ) by C(v, h) = g where £ 6 V is 
the unique solution to 

<a(|, xjr) + b(v, z\ xj/) + h(z*. u, if) - kH (h, x)r, =0 for xf = (0, *) € V. 

Then, since H l (£2) is embedded compactly to L 4 (fi), Lemma 1 implies that C is 
compact. Thus, by the Banach closed range theorem and the Riesz-Schauder theo- 
rem, e'(jc*)(u, h) is sutjective if and only if ker(G*) — {0} [4], which is equivalent 
to (21). □ 

4. Augmented Lagrangian method 

In this section we discuss applications of the augmented Lagrangian method for 
the constrained minimization problem (17). The augmented Lagrangian method 
[8], [15] is based on an equivalent formulation of (17): 

c , 

minimize J(x) + - \e(x)\ Y over x e X and g e C . (24) 

subject to e(x) = 0, where x = ((«, T - T 0 ),g) 6 X = V x L 2 {T\) and c > 
0 is the penalty parameter. The augmented Lagrangian algorithm [8], [15] is the 
multiplier method applied to (17), i.e., it involves a sequence of minimizations of 
the functional 


L Ct (r. A.*) = J(x) + ( A* , e(x)) + c - \e(x)\ 2 Y 
subject to g e C , 


(25) 
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(25) 


where the multiplier sequence (A.*} in '¥' is generated by the first order update 

X k+ '-X k -r c k e(x k ). (26) 

for k > 1. where .v * is a minimizer of L Ci (-. X k ) and we assumed that Y' = Y . 
otherwise each element in K* has its Riesz representation. To carry out this it- 
eration. a sequence of monotonically nondecreasing, positive real numbers { cy- 1 . 
c\ > Co > 0 and a start up value A 1 for the Lagrange multiplier for the equality 
constraint e(x) = 0 need to be chosen. The convergeneefesults of the augmented 
Lagrangian method for the infinite dimensional optimization problem are estab- 
lished, for example, in [9], [16]. The augmented Lagrangian method is a hybrid of 
the penalty method (i.e., X k = 0) and the Lagrange multiplier method (i.e.. c* = 0) 
and combines good properties of both. It overcomes the difficulty of the penalty 
method which requires a large value of c*. The augmented functional L lt (x. A 4 ) 
is locally strictly convex provided that A* is sufficiently close to A* and the second 
order optimality condition 

Ll{x\ A*)((u. h). ( v . h)) > a (|u| 2 v + \h\l) 

for all (u, h) e X satisfying e'(.r')(ti. h) = 0 . 

for some a > 0, is satisfied. Here, Lq (x*. A*) denotes the bilinear form that char- 
acterizes the second derivative of Lo(x, A) = 7(.r)-|-(A. e(x)) with respect to .r at 
x*. That is, it is not necessary that the cost functional J be (locally) convex, which 
is required for convergence of the multiplier method. The algorithm (25)-(26) has 
been successfully applied to parameter estimation problems in elliptic PDEs [10], 
[12] and optimal control problems for 2-D incompressible Navier-Stokes [4], The 
first order update (26) provides Q-linear convergence of the iterates .r* in X. In [ 1 1 ] 
we have investigated a second order update scheme for the augmented Lagrangian 
method. In what follows we assume that g* eint (O- Thus. (19) reduces to t?'(.r*) 
is surjective (see Lemma 6). Hence the necessary condition implies that 


L' c (x*, A*) = 0 and. e{x M ) = 0, 


(28) 


for all c > 0. An algorithm proposed in [11] applies Newton’s method to (28). 
The resulting algorithm is stated as: given a current iterate (.r. A) the next iterate 
(,t+, A+) satisfies 


( K(x, A) e'{x)* \ / x + -x 
V e\x) 0 ) V A + - A 


L' c (x. A) 
e(x) 


(29) 


Note that 


and 

L"(x, A) 


L' c (. t,A) = L' 0 (x, X + c e(x)) 

Lq (x, A + c e(x)) + c {e'(x)( ), e'(x)(-)) . 


(30) 
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Consequently, suppose |(x, A) — (x*, A*)| is sufficiently small. Then it follows 
from (27) [11] that L"(x, A) is coercive on A x Y . Thus equation (29) can be re- 
garded as a general Stokes equation. Following an argument due to Bertsekes [2], 
we can avoid forming L" during the iteration. From the second equation of (29) 
we have e'(x)(x+ — x) = —e(x). Thus the first equation can be written as 


Lq(x, k + ce(x))(x+ - x) + (?'(x)’(A + - (k + ce(x))) = -L' 0 (x, k +ce(x)) 
and hence (29) is equivalent to of 


L"(x,k) e'(x)* 

e\x) 0 




(31) 


where A. = k -I- c e(x). 

Note that A is nothing but the first order update of the Lagrange multiplier if the cur- 
rent iterate x minimizes L c (x, k). Equation (31) is more advantageous than (29) 
since the squaring term ce'(x)*e'(x) is absorbed and less calculation is involved. 
If we define a matrix operator S on X x Y by 


Remark. A 
duced to the 
reduction o 
ample, if jc 
A ssume tha 
ciently sma 
Q-quadratic 

5. Hybrid 

In this sectic 
retical result 
dimensional 

minimi/ 


subject to (2 
interior of f 


Six, k) 


L"(x,A) e\xY 
e'(x) 0 


then it follows from (27) that S(x*, A*) is boundedly invertible. Thus, if (x, A) is 
sufficiently close to (x*, A*), then equation (31) has a unique solution. We sum- 
marize our discussions as 


Algorithm 1 

(1) Choose A 1 e Y, c > c > 0, and set c — c — c, k = 1. 

(2) Determine x = ((u,T — T 0 ), g) e V x C such that 

L c (x,k k ) < L c (x*,k k ) = fix'). 

(3) Set A = A* + c e(x). 

(4) Solve for (x + , A + ) € X x Y\ 


S(x , A) 


(5) Setx^+i = x+ and A* +l 
set k = k + 1 and go to (2). 


*+“•* \ / L’ 0 ix,k) 

k+ - A / \ e(x) 

= A + . If the convergence criterion is not satisfied then 
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Remark. A variant of Algorithm l is obtained by skipping step (2). Then it is re- 
duced to the Newton method applied to equation (28). Step (2) implies a sufficient 
reduction of the merit functional (the augmented Lagrange functional). For ex- 
ample, if x — (z. g) minimizes L c (-, A.*) over V x C then step (2) is completed. 
Assume that (19) and (27) hold. It is proved in [11] that if | A. 1 — A. “ | >- is suffi- 
ciently small, then Algorithm 1 is well-posed and (.r*. X k ) converges to (,r ’ . X') 
Q-quadratically. 

5. Hybrid method and test example * 

In this section we present an example to demonstrate the applicability of our theo- 
retical results and proposed algorithm. We consider the optimal control of the two 
dimensional stationary thermally driven cavity flow 

minimize J{g) = \ f Q |ti - u d \ 2 + |T - T d \ 2 dx -F § |g - Toi^ ;(ri) 
over g e L 2 (T i) 

subject to (2), where £2 = (0, L) 2 , Ti = (,t| , 1). 0 < x\ < L and To is the relative 
interior of f — IY On Tq we have the Dirichlet boundary condition: 

7(0, x 2 ) = T(L, x 2 ) = T 0 , 0 < x 2 < L 

T(.t,,0) = T 0 -50 min(l, = (2.5 - j{x, -2.5|), 0 < x, < L . 

with a reference temperature 7 0 = 1350° AT, and the temperature on the substrate 
(xi , 0), -yL < .t| < yL is 1300° AT The desired state (u d , Tu) appearing in 
(32) is chosen as follows, uj = L u° where (u°, T°) is the solution pair for the 
flow with L = 1 and g — T 0 and Tj — T' where («’, 7 1 ) is the solution pair 
for the flow with L = 5 and g = T 0 . Here, our numerical calculation strongly 
indicates that the solution to (2) is unique and (dynamically) stable for the range 
of physical parameters that we chosen for our calculation. 

The problem is scaled so that the velocity field u has dimension cm I sec and 7 
has dimension K° . The constants v, k are chosen for P 2 at 3 atm pressure and are 
given by v = . 155 and k = .1 10. H, /3 and L are set to be H = 100, fi — = 

5.5 x 10~ 5 and L = 5 respectively. Figures 1 and 2 show the vector field of «° 
and u 1 , respectively. It can be observed that u° has more vertical transport than 
u' does and u 1 is confined in the two bottom corners. Hence the cost-functional 
(24) is formulated so that the thermal control g on the top T i increases the vertical 
transport of flow with L = 5 while retaining the temperature distribution 7 1 at the 
reference temperature g = Tq. 
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Figure 1. The desired flow. 



Figure 2. The uncontrolled flow. 

The second order augmented Lagrangian method (Algorithm 1) described in 
Section 4 is used to solve problem (32). To obtain a good starting value for 
the Lagrange multiplier in Algorithm 1, we employed a few steps of the gradi- 
ent method. The gradient of the cost functional can be calculated by the adjoint 
equation as follows. Consider the cost functional J (z, g ) subject to E(z, g) = 0 
where z € V and g e U. Assume that V and U are Hilbert spaces, J(z, g) : X = 
VkU -> R is Frechet differentiable and E(z, g ) : X — ► Y is continuously Frechet 
differentiable in a neighborhood of xq = (zo, go)- Suppose that E z {xq) e £(V, Y), 
the F— derivative of E with respect to z at xq, has a bounded inverse. Then by the 
implicit function theorem, there exists a unique C 1 mapping ^ defined in a neigh- 
borhood N of go in U, V : U -*■ V such that 4>(go) = Zo and £(fl/(g), g) = 0 
for g € N. Then the F — derivative d of /('l'(g), g) with respect to g in N exists 
and is given by 

d = J g + E' g f (34) 


wh 


We 
als 
me: 
for 
[3]. 
•n [ 
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where \j/ e Y satisfies the adjoint equation 


E:# + J z =0. (35) 

Here we assume that V, U and Y are identified with their dual spaces. In fact, if 
u = H/'(/i), h 6 U then 


and u e V satisfies 


Hence, for h e U 


(d. h)u = ( J ,. , h)u + (J., u)v 


E-u + E.,h — 0 . 


(d, h) u = (J x , h) u - v) v = (•/,., h)u + (X. E-h) v = (J g + ETlt. h) 

which implies (34). The projected gradient method can be written as: 

Algorithm 2 

(1) Choose °\ e U and set k = 1. 

(2) Let ~k be a solution of E(z, gk) = 0, t fr k be the solution of E : (z k , gkV'fa + 
Jzkk, 3k) = 0 and set d k = J s {z k , g k ) + E s (z k , gkY^k- 

(3) Set dk = 7 s (zt,.ft)and determine a* > Osuch that 7(^(g a ), g a ) is minimized 
where g a = Projcigk ~ a d k ). 

(4) Set g k+ 1 = g k — a k d k . If the convergence criterion is not satisfied, then set 
k — k + 1 and go to Step (2). 

In our specific example equation (35) is written as 

-v Aa — u ■ VA. + X'Vu - (T — T 0 ) + Vq + u — Uj = 0 

V ■ X = 0 and X|p - 0 

° (36) 

—k A fi — « • V/z — — Xt + T — Tj = 0 
To 

yz = 0onfo and n ■ V q. + xH fx = 0 on F| , 

where >/r = (a, /j.) e V = V 0 x V { , q e L : (Q) and 

£* i/r = — /erf /z on Ti . (37) 

We used the mixed-finite element method [7] based on the Legendre polynomi- 
als to approximate problem (l)-(2) numerically. Detailed discussions about the 
method are given in [13], In our implementation we calculated the adjoint system 
for the approximated problem and solved equations (2) and (36) using GMRES 
[3], The specific implementation of GMRES applied to our example is described 
in [13] in detail. Concerning the divergence-free constraint V u = 0, we employed 
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the feasible method, projecting the first equation onto the divergence-free space V 0 
as in [4], [6], [13]. The line search in Step 3 of Algorithm 2 was performed by the 
linearization of the constraint £(z, g) = 0 at ( z k , gk ) since the cost functional is 
quadratic. That is, if v k is the solution to E z (z k , gk)v k -F E s {z k , g k )d k = 0 then 
a > 0 is chosen so that J(z k — a v k , g k — ct d k ) is minimized. 

For this specific example, three steps of Algorithm 2 were performed. We then 
set x 1 and A. 1 as x 1 = (24, g 4 ) and A. 1 = ^3 for Algorithm 1 without Step 2. The 
matrix operator S(x, A) was calculated for the approximated system and the result- 
ing linear equation (31) was again solved by GMRES. ^ 

The calculations were performed using a 20 x 20 Cartesian product of Legendre 
polynomials, choosing c = 1 and g\ = 0. Algorithm 1 was terminated after three 
iterates, since the necessary and sufficient optimality condition (28) was satisfied 
within a residual norm of 1 x 10 -7 . 

We may compare this rapid convergence of the hybrid method with the results 
of using either algorithm by itself. Algorithm 2 did not fully converge after 50 
iterates. Algorithm 1, with the start-up x 1 = ((«', T 1 ), 0) 6 X and A 1 = 0, 
also failed to converge. Thus, the use of the hybrid method combining the gradient 
method and the second-order augmented Lagrangian was essential for the success 
of our numerical calculations. Figure 3 shows the iterates g k (the first three curves 
from top to middle) for Algorithm 1 and the (calculated) optimal control g * (the 
lowest curve). The iterates for Algorithm 1 are not shown because they coincide 
with g * within the accuracy of the plotting. Figure 4 shows the resulting vector 
field «* which corresponds to g”. It shows clearly that the vertical transport of 
flow is increased. 



Figure 3. Optimal control iterates. 
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Figure 4. The controlled flow. 
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